Interaction of streamers in air and other oxygen— nitrogen mixtures 
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The interaction of streamers in nitrogen-oxygen mixtures such as air is studied. First, an efficient 
method for fully three-dimensional streamer simulations in multiprocessor machines is introduced. 
With its help, we find two competing mechanisms how two adjacent streamers can interact: through 
electrostatic repulsion and through attraction due to nonlocal photo-ionization. The non-intuitive 
effects of pressure and of the nitrogen-oxygen ratio are discussed. As photo-ionization is experimen- 
tally difficult to access, we finally suggest to measure it indirectly through streamer interactions. 
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Streamer discharges are fundamental building blocks 
of sparks and lightning in any ionizable matter; they are 
thin plasma channels that penetrate nonconducting me- 
dia suddenly exposed to an intense electric field. They 
propagate by enhancing the electric field at their tip to a 
level that facilitates an ionization reaction by electron im- 
jfl . Streamers are also the mechanism underlying 
I) 0, ; these are large atmospheric discharges 
above thunderclouds that, despite being tens of kilome- 
ters wide and intensely luminous, were not reported until 
1990 [6| . Although the investigation of streamers concen- 
trates mainly in gaseous media, they have also been stud- 
ied in dense matter, such as semiconductors Q and oil 0. 
Streamers have also received attention in the context of 
Laplacian-driven growth dynamics Q and a strong anal- 
ogy with viscous fingering, in particular Hele-Shaw flows 
UM, has been established. 



Both in laboratory 11| and in nature [12j |, streamers 
appear frequently in trees or bundles. As their heads 
carry a substantial net electrical charge of equal polar- 
ity that creates the local field enhancement, they clearly 
must repel each other electrostatically which probably 
causes the "carrot"-like conical shape of sprites. On 
the other hand, recent sprite observations tl3| as well as 



streamer experiments (|l4|. Fig. 7, [ill ]. Fig. 6) also show 
the opposite: streamers attract each other and coalesce. 

Up to now, streamer interactions have not been stud- 
ied much theoretically, and streamer attraction has not 
been predicted at all. In coarse grained phenomenologi- 
cal models for a streamer tree as a whole [151 ] , the repul- 
sive electrostatic interaction between streamers is taken 
into account. In a more microscopic, but still largely 
simplified model, Naidis studied the corrections to 
the streamer velocity due to electrostatic interaction with 
neighboring streamers. In [l7j], two authors of the present 
letter have studied a microscopic "fluid" model for a pe- 
riodic array of negative streamers in two spatial dimen- 
sions, where they show that shape, velocity and electro- 
dynamics of an array of streamers substantially differs 
from those of single streamers due to their electrostatic 
interaction, but attraction or repulsion were excluded by 
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FIG. 1: Two negative streamers in nitrogen at atmospheric 
pressure advancing downwards and repelling each other; 
shown are surfaces of constant electron density in an advanced 
state of evolution within a constant background field. See the 
text for details. 



the approach. Due to the difficulty to represent this mul- 
tiscale process Q in a numerically efficient manner, only 
recently it has become possible to simulate streamers in 
full 3D [3 We here present a numerical method to 
handle this problem, and we apply it to the interaction 



of streamers in complex gases like air where a nonlocal 
photon mediated ionization reaction has to be taken into 
account. We find that when varying gas composition and 
pressure, streamers can either repel or attract each other. 
The transition occurs in an unexpected manner, and is 
not simply determined by an ionization length. 

The photon mediated ionization reaction in air and 
other nitrogen oxygen mixtures is experimentally not eas- 
ily accessible, but forms a basic ing redient of the present 
theory of streamers in air [3 M, HE rH ; all these 
simulations are based on the single experimental mca- 
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surement of Penney and Hummert in 1970 2J|. Our theo- 
retical results suggest that this reaction could be deduced 
from experiments on coalescence or repulsion of adjacent 
streamers as a function of pressure and gas composition. 

We study the minimal streamer model 25j extended by 
the nonlocal photo-ionization reaction characteristic for 
nitrogen oxygen mixtures like air. It consists of continu- 
ity equations for electron and ion densities n e ,+ coupled 
to the electrical field E that is determined by the poten- 
tial on the outer boundaries and space charge effects 

d t n e = V • (rieMeE) + D e V 2 n e + 5, ; + S ph , (1) 
d f n + = Si + S ph , (2) 
e V-E = e(n+-n_), E = -V</>. (3) 

Here \x e is the electron mobility, D e is the electron diffu- 
sion coefficient and e is the elementary charge. Ion mobil- 
ity, much smaller than electron mobility, is neglected. To 
fully focus on the influence of photo-ionization, electron 
attachment on oxygen is here neglected as well, and we 
use transport parameters for pure nitrogen as in previous 
work 0, [26} . The source terms for additional electron- ion 
pairs are the local impact ionization Si in Townsend ap- 
proximation, Si = n e ^ e |E|a(|E|) = n e ^ e |E|a e _£ ' t " / ' E ', 
where ao is the ionization cross section and Eq is the 
threshold field, and the nonlocal photo-ionization accord- 
ing to the model for oxygen- nitrogen mixtures developed 
by Zhclezniak et al. |27j 
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with A(p) — p q / (p + p q ). Here it is assumed that accel- 
erated electrons excite the b 1 ^, b a H+ and C4S+ states 
of nitrogen by impact with a rate £Si where Si is the lo- 
cal impact ionization rate and £ a proportionality factor. 
These nitrogen states can deexcite under emission of a 
photon in the wavelength range 980 — 1025 A that can 
ionize oxygen molecules [3, [27j . The absorption length 
of these photons by oxygen is obviously inversely propor- 
tional to the oxygen partial pressure po 2 ■ Introducing the 
oxygen concentration in the gas as rj = po 2 /P> the absorp- 
tion function of photoionizing radiation h(pr) is charac- 
terized by the two length scales pr m i n 1=3 380 /im • bar/7/ 
and pr ma x ~ 6.6 /im • bar/r/. Above a critic al g as pres- 
sure that we take as p q = 60 Torr = 80 mbar [22j , the ex- 
cited nitrogen states can be quenched by collisions with 
neutrals, hence suppressing photo-emission; this is taken 
into account through the prefactor A(p). 

A major difficulty in the theoretical study of streamer- 
streamer interactions is the development of an efficient 
numerical code for streamer simulations with their inner 
multiscalc structure in three dimensions. Here we present 
a locally adaptive and parallelizable approach to this 
problem. Our method has in common with the one de- 
scribed in [l8j that cylindrical coordinates and a uniform 



grid in the angular dimension are used, but in the r, z- 
projection we apply the grid refinement scheme of [26| to 
resolve better the thin, pancake-like shape of the space- 
charge layer. Furthermore, to allow the parallel solution 
of the Poisson equation in multi-processor machines we 
perform a Fourier transformation in the angular coordi- 
nate 9, (f> k {r 7 z) = J2n=o <M r > z i n )e~ lke ™ , where N is the 
number of grid cells in the 8 direction and 9 n = 2im/N. 
For each separate mode k, a two-dimensional Hclmholtz 
equation has to be solved: 
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where a tilde~represents the Fourier transform of a quan- 
tity and \wk\ 2 = — cosfcA#). For each Fourier 
mode, we apply the refinement algorithm described in 
[111 , which is trivially generalized to solve the Helmholtz 
instead of the Poisson equation. The advantage of solv- 
ing the electrostatic part of the problem in the Fourier 
domain is that each of the Fourier modes is decoupled 
from the rest and hence it can be solved in parallel in 
a multi-processor computer. Once the electric field is 
calculated in Fourier space, the field in real space is de- 
rived through an inverse Fast Fourier Transform; and the 
convection-diffusion-reaction system {l])-© is integrated 
as detailed in [26J. The photo-ionization term is com- 
puted in a Helmholtz PDE approach as described in [13] , 
and in the angular direction with a scheme of Fourier 
transformations and parallel solving that is completely 
equivalent to the one applied to the Poisson equation. 

We now use our model and numerical algorithm to 
study the interaction between two streamers. We fo- 
cus on the influence of (i) the pressure in air and of 
(ii) the oxygen-nitrogen ratio at standard pressure; stan- 
dard temperature is always assumed, (i) The compari- 
son of streamers at different pressures p relies on scaling 
lengths, times etc. with appropriate powers of p. These 
similarity laws are strictly valid in the minimal streamer 
model 0, ; they are broken by photo-ionization (g]) 
when the pressure reaches the quenching pressure p q ; for 
P ^ Pqi photo-ionization at unchang ed g as composition 
is increasingly suppressed like A(p) [20|. Similarit y o r 
Townsend scaling is further discussed in Sect. 1.2 of [29j . 
(ii) The 02:N2 ratio changes the local photoionization 
rate which is proportional to it, while the absorption 
lengths are inversely proportional to it. 

First, negative streamers are investigated, because 
they propagate even in the absence of photo-ionization. 
Then it is shown that positive streamers behave sim- 
ilarly. In all simulations, we use two identical Gaus- 
sian seeds of reduced width (according to Townsend scal- 
ing) p ■ w — 73.6 /im • bar and amplitude p~ 2 • n e ^ max 
1 .4 • 10~ 2 fj,m~ 3 ■ bar -2 separated by a reduced distance 
of p ■ d = 230/im • bar as an initial condition. They are 
exposed to a homogeneous and constant background elec- 
tric field Eb/p = 80kV/(cmbar) in the positive z direc- 
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FIG. 2: Evolution of the space charge layers of two adja- 
cent negative streamers at different pressures in air. (a) 
p — 0.07mbar, (b) p = lbar, (c) p = 50 bar, (d) p — * oo. (a) 
corresponds to A(p) w 1 and (d) to A(p) = 0, which is equiv- 
alent to pure nitrogen since there is no photo-ionization (even 
though in pure nitrogen the length-scales formally diverge). 
The axes show reduced lengths p-x, p-z to exhibit similarities. 
In the upper right corners, the two photo-ionization lengths 
are inserted as vertical bars. A multiplicative factor is used 
where the lengths are too long to fit into the figure. 

tion. The streamers on reduced length, time and density 
scales are similar if photo-ionization is neglected. 

The simulations were performed with an angular reso- 
lution of A6 = 27r/64 (N = 64), but we also checked that 
the results are stable when we double the number of an- 
gular grid cells. Figure 1 shows a surface of equal electron 
density for the discharge in nitrogen at normal pressure 
and temperature at time t — 1.56 ns. Figures 2-4 show 
the space charge layers at the streamer heads (more pre- 
cisely, the half maximum line at the respective time) in 
the plane intersecting the two streamer axes in timesteps 
of 0.12 ns • bar; the first snapshot is at time 0.84 ns • bar. 

Figure 2 demonstrates the effect of pressure change 
for similar streamers in artificial air, i.e., in an oxy- 
gen nitrogen mixture of ratio 20:80, corresponding to 
rj = 0.2. The panels show (a) air at a pressure of 
0.05 torr w 0.07 mbar, which is the atmospheric pres- 
sure at approximately 70 km height, where sprites are 
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FIG. 3: Evolution of the space charge layers of two adjacent 
streamers at atmospheric pressure, but different concentra- 
tions of 2 : T) = po 2 /p = 10~\ 10 -2 , 10 -3 , 1(T 4 in panels 
(a)-(d). The two photo-ionization lengths are again inserted 
as vertical bars in the upper right corners. 

frequently observed, (b) air at 1 bar, (c) air at 50 bar, 
and (d) pure nitrogen which in our model is equivalent 
to air at infinite pressure. 

The first observation is that two streamers do not al- 
ways repel each other. Rather there are two qualitatively 
different regimes: at a pressure of 50 bar and above, when 
the fast collisional quenching of N 2 molecules suppresses 
photoionization, the streamers repel each other electro- 
statically due to the net negative charge in their heads. 
However, at atmospheric pressure and below, the stream- 
ers attract each other: a cloud of electrons is created 
between the two streamers which eventually makes them 
coalesce into a single, wider one. As it has been suggested 
that the photo-ionization length could determine length 
scales in the streamer head [30j, the two reduced photo- 
ionization lengths are indicated by the vertical bars in the 
upper right corners, they are the same in all panels. Ob- 
viously, there is no simple relation between these lengths 
and the observed repulsion or attraction; the interaction 
is rather governed by the quenching prefactor A(jp). 

Figure 3 demonstrates the effect of the oxygen con- 
centration. For fixed atmospheric pressure, the rel- 
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FIG. 4: Evolution of the space charge layers of two adjacent 
positive streamers in air at atmospheric pressure at 1 bar (a) 
and 50 bar (b). 



ative oxygen concentration rj = po 2 /p is reduced to 
1CT 1 , 10" 2 , 10" 3 , 10" 4 in panels (a)-(d). As the photo- 
ionization lengths scale inversely with the oxygen con- 
centration, these lengths increase by factors 10 from one 
panel to the next while the prefactor decreases with a 
factor 10 _1 . The attraction between the streamers de- 
creases with decreasing oxygen concentration until it is 
not visible anymore in panel (d) for 77 = 1CP 4 . 

Figure [4] shows that the interaction between two posi- 
tive streamers is qualitatively similar to that of two neg- 
ative streamers, since the same two competing phenom- 
ena are present. Initial and boundary conditions are the 
same, and panels (a) and (b) show positive streamers in 
air at atmospheric pressure and at 50 bar. The positive 
streamers take longer to start but after a relatively short 
time they merge. 

Due to the nonlinear and nonlocal nature of streamer 
interactions and to the many dimensions of our param- 
eter space, it is difficult to develop a simple prediction 
for streamer merging. Nevertheless, the following re- 
marks can be made. First, to produce enough ioniza- 
tion between the streamers, at least one of the absorp- 
tion lengths of photo-ionization must be larger than or 
comparable to the streamer distance. Second, this photo- 
ionization is only amplified to a level comparable to that 
on the streamer head if the field is enhanced between 
the streamers; hence the streamer distance should not 
be much larger than the streamer radius. Finally, an ap- 
proximation for the relative photo-ionization level can be 
extracted from (j4|) as the maximum of the instantaneous 
rate of photo-ionization on the middle axis produced per 
impact ionization event in the tip of a given streamer, 
namely (3 — ^A(p)h(pd/2)/Trp 2 d 2 . Preliminary results 
show that this parameter influences the electron densi- 
ties in the axis at early stages of the evolution. We have 
not found a deterministic law but coalescence is favoured 
for large f3 and occurs always if j3 > 10 -11 ^m -3 • bar -3 . 

We have developed a code to study the interaction of 
streamers in full 3D space, and we have studied the ba- 



sic processes that govern the interaction of two adjacent 
streamers. To focus on the underlying physics of photo- 
ionization we have neglected electron attachment and 
nontrivial electrode geometries. Further steps are needed 
to successfully predict the outcome of experiments and 
observations; we mention needle electrodes and the non- 
alignment of the streamer heads. But our results show 
that for a given pressure p, electric field Eq, oxygen- 
nitrogen ratio 77 and initial seeds, there is a threshold 
distance d* below which two streamers coalesce. This 
distance, which is experimentally accessible, would be 
an indirect measure of the frequency of photoionizing 
events. Hence we believe that fully three-dimensional 
calculations of streamer dynamics will provide a suitable 
test-case for streamer models as they predict easily ob- 
servable behavior like attraction, repulsion or branching. 
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